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FOREWORD 


This research was supported under contract NAS1-10277 with the 
National Aeronautics and Space Administration’s Langley Research Center. 
The technical monitor on the contract was H. Harris Hamilton. Previous 
results, without the effect of variable entropy at boundary layer edge, 
were published in NASA CR 111921 and NASA CR 111922 in August, 1971. 

This report gives a description of the theoretical approach and changes 
for the computer program to include the effects of variable entropy at 
boundary layer edge. 
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CALCULATION OF HEAT TRANSFER ON SHUTTLE TYPE CONFIGURATIONS 


INCLUDING THE EFFECTS OF VARIABLE ENTROPY AT BOUNDARY LAYER EDGE 

by Fred R. DeJarnette 
North Carolina State University 

October, 1972 
SUMMARY 

A relatively simple method is presented for including the effect of 
variable entropy at the boundary-layer edge in a heat transfer method devel- 
oped previously. For each inviscid surface streamline an approximate shock- 
wave shape is calculated using a modified form of Maslen's method for invis- 
cid axisymmetric flows. The entropy for the streamline at the edge of the 
boundary layer is determined by equating the mass flux through the shock 
wave to that inside the boundary layer. Approximations used in this tech- 
nique allow the heating rates along each inviscid surface streamline to. be 
calculated independent of the other streamlines. 

The shock standoff distances computed by the present method are found 
to compare well with those computed by Maslen's asymmetric method. Heating 
rates are presented for blunted circular and elliptical cones and a typical 
space shuttle orbiter at angles of attack. Variable entropy effects are 
found to increase heating rates downstream of the nose significantly higher 
than those computed using normal-^shock entropy, and turbulent heating rates 
increased more than laminar rates- Effects of Reynolds number and angles 
of attack are also shown. 



INTRODUCTION 


A relatively simple method was developed in Reference 1 to calculate 
laminar, transitional and turbulent heating rates on arbitrary blunt-nosed 
three-dimensional bodies at angle of attack in hypersonic flow, Inviscid 
surface streamlines were calculated from Euler's equation using a prescribed 
pressure distribution,, Heating rates were determined along a streamline by 
applying the axisymmetric analog (small cross-flow approximation) for three- 
dimensional boundary layers to solutions of the axisymmetric boundary-layer 
equations. In this approximation the distance along a streamline is inter- 
preted as the distance along an equivalent axisymmetric body, and the scale 
factor (or metric coefficient) for the surface coordinate normal to the 
streamline (which is a measure of the streamline divergence) is interpreted 
as the radius of the equivalent axisymmetric body. Thus, each inviscid sur- 
face streamline corresponds to a different equivalent body of revolution at 
zero incidence. In Reference 2 this method was shown to yield accurate 
laminar heating rates on blunted circular and elliptical cones, and a typi- 
cal delta-wing space shuttle orbiter at angle of attack. 

In the method discussed above the; properties at the edge of the bounda- 
ry layer were determined from the surface pressure and the entropy aft of a 
normal shock wave (typical blunt body assumption) for an ideal gas or equi- 
librium air. While this assumption is valid in the nose region of a very 
blunt bodies, its accuracy decreases as the ratio of distance from the stag- 
nation point to nose radius increases. As the boundary layer along the sur- 
face grows, more and more mass is entrained into the boundary layer and it 
is therefore possible for the boundary layer to swallow the entropy layer. 
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This means that the streamlines which passed through the nearly normal por- 
tion of the bow shock wave are now inside the boundary layer. Although the 
pressure at the edge of the boundary layer is still nearly the same as that 
on the surface, the entropy of the streamline at the edge of the boundary 
layer can be quite different from the normal-shock entropy. As a result, 
the corresponding heating rate is generally higher than that calculated us- 
ing the normal-shock entropy. 

This report develops an approximate method to include variable entropy 
at the edge of the boundary layer in the method of References 1 and 2 for 
calculating heating rates. Modifications needed to combine this addition 
with the computer program given in Reference 3 are also described. Results 
are presented for hypersonic flow over blunted circular and elliptical cones 
and a typical space-shuttle orbiter at angle of attack. 


SYMBOLS 


*1 * 2 * **3 


B 

B 

C 


coefficients for geometry of space shuttle orbiter, see eq. 
(76) 

ratio of body principal radii of curvature, R^ ^ 

ratio of shock principal radii of curvature, R^/R-^ 
parameter used in eq. (18) 


C' 


D 


A A 


V e 0’ e n 


dC/d& 

A 

constant defined by eq. (59) 

unit vectors in shock-oriented coordinate system, eqs. (2), 
(10) , and (11) 
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e- 

x 


» 


e- 

y 



E 

f 

G 2 

h 

h 

h 

s 

h 

c 

h E 

h w 

h 2 

Ah 

H 

H 

s 

L 

M 

ri 

ri 


P 

Ap 

% 

q 

r 

Re/ft 


unit vectors in wind-oriented Cartesian coordinate system 
constant defined by eq. (60) 
body radius, ft, 

1 - sin 2 r 

scale factor in 3-direction on body, ft. 

scale factor in 3-direction on shock, ft. 

scale factor in £ direction on shock, dS = h d£, ft. 

s 

2 

heat transfer coefficient, Btu/ft -sec-°R 

8 2 2 

.reference enthalpy, 2,119 x 10 ft /sec • 

2 2 

wall enthalpy, ft /sec 

2 2 

enthalpy aft of normal shock, ft /sec 

2 2 

enthalpy parameter defined by eq. (46), ft /sec 
form factor 6*/8 

m 

2 2 

stagnation enthalpy, ft /sec 
body length, ft. 

Mach number 

distance normal to a streamline, ft, 

straight line coordinate normal to shock wave and toward body 
ft, 

2 

pressure, lb/ft 

2 

pressure parameter defined by eq, (45), Ib/ft 

2 

heat-transfer rate at wall, Btu/ft -sec. 

distance norma! to 3 ® constant lines on shock, dq * h d3, ft 
radius, ft, 

frees tream Reynolds number per foot 
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*N 


R T,b’ R H,b 

-> 

R 

R 


freestream Reynolds number based on nose radius 
nose radius, ft. 

body principal radii of curvature at stagnation point, ft. 
position vector, ft. 

shock radius of curvature along I = constant line, see eq. 
(43), ft. 


s 

s 

T' 

u,v,w 


V 

* 

W f 

x,y,z 


x.y.z 


X 


* 


3 

3 

Y 

r 


shock principal radii of curvature at stagnation line, ft. 
distance along inviscid surface streamline, ft. 
distance along 3 = constant line on shock, ft. 
temperature, °R 

velocity components in £■> n, and 3 directions, respectively, 
ft/sec 

velocity magnitude, ft/sec 
velocity vector, ft/sec 
weighting function 

body-oriented coordinate system, see Fig. 14 
Cartesian coordinates in wind-oriented coordinate system with 
origin at the stagnation line of shock wave (see Fig, 1) 
distance defined by eq. (77), ft. 

distance from nose to beginning of body segment, ft. 
body coordinate normal to inviscid surface streamline 
shock coordinate normal to \ lines 
ratio of specific heats 
shock-wave angle, see eq . (3) 
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m 

* 

tr 

w 


A 

l 


P 

a 


<f> 


4 > 

'V 


body angle relative to V^, see eq. (B13) 
boundary-layer thickness, ft, 
boundary-layer displacement thickness, ft, 
transformed displacement thickness, ft, 
shock standoff distance at stagnation line, ft. 
parameter defined by eq. (25), ft, 
momentum thickness, ft, 
transformed. momentum thickness, ft, 
wall enthalpy ratio, h w /H g 
- 

exponent defined by eq, (54) 

coordinate along shock surface 

3 

mass density, slug/ ft 

angle measured on shock wave, see eq. (4) and Fig. 2 

body circumferential angle with <j> = 0 in windward plane of 

symmetry 

shock circumferential angle, ? “ tan ^ (z/y) 

second stream function, see eq. (7) 

first stream function, see eq. (7), slugs/sec 


Subscripts: 


b 

B.L. 

e 

lam 


body 

boundary layer 

edge of boundary layer 

laminar 
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s 

sh 

turb 

w 

e 

2 


stagnation point 
at shock wave 
turbulent 
wall 

evaluated at x- = e / (211^) 
aft of normal shock 
freestream conditions 
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ANALYSIS 


Any highly accurate method which accounts for variable entropy at the 
edge of the boundary layer would require a complete solution to the three- 
dimensional inviscid flow field. Such a technique, however, would destroy 
the simplicity and short computational time that are characteristics of the 
computer program of Reference 3. Therefore, an approximate, but reasonably 
accurate, variable entropy method is developed here which is consistent with 
the other approximations involved in Reference 2. 

In order to determine the entropy at the edge of the boundary layer, the 
position where the edge streamline passed through the shock wave must be 
determined. Then the slope of the shock wave at this position is all that 
is needed to calculate the entropy of that streamline. Therefore, a complete 
knowledge of the inviscid flow field is not needed. 

For axisymmetric inviscid flows, Maslen (Ref. 4) developed a simple in- 
verse method based on a von Mises transformation coupled with a simple 
approximate integral of the lateral momentum equation. Results for several 
examples were found to compare favorably with those of more exact methods. 

In particular, the shape of the shock wave and surface pressure distribution 
(except near the stagnation point) were shown to be very accurate. This 
method was extended to asymmetric flows in Reference 5; however, only re- 
sults for the stagnation region of blunt bodies and conical flows were 
presented in that report. The technique discussed therein for solving the 
general three-dimensional inviscid flow field is much more complicated than 
the method developed for axisymmetric flows (Ref. 4). The method described 
below represents an extension of the axisymmetric Maslen method to general 
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three-dimensional flows; however, it is much simplier than Haslen's asymmet- 
ric method (Ref. 5). The stagnation region will be considered first, since 
the results can be compared with Maslen’s numerical results. 

Stagnation Region 

As noted by Maslen (Ref. 4), the approximations employed in his method 
are not strictly valid in the stagnation region; however, for axisymmetric 
flows it was found to give a reasonably accurate value for the shock stand- 
off distance. For asymmetric flows the requirement is more severe because 
in addition to a reasonable standoff distance, results should be independent 
of the direction from which the stagnation line is approached. Maslen’s 
asymmetric solution (Ref. 5) yielded different standoff distances for dif- 
ferent directions of approaching the stagnation line. The solution obtained 
below yields a unique value for the shock standoff distance. 

For the stagnation region, it is convenient to use wind-oriented 
Cartesian coordinates x, y, z whose origin is located at the normal position 
on the shock wave with x in the direction of (see Fig. 1). Following 

Maslen, the shape of the shock wave near the stagnation line is represented 
by an elliptic paraboloid 



2Rn 2^ 


where the x, y plane is a plane of symmetry. This equation represents the 
shock wave by a portion of an ellipsoid for x small since the equation is 
an ellipse in y and z for x = constant. At the origin, the shock radius 
of curvature is Rji in the x, y plane and R^ in the x, z plane. 
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The unit (outer) normal to the shock surface is illustrated in Fig. 2 


and is given by the relation 


e_ 

n 


A 


- e. 


x 


+ 


-Z-a. + — e- 

Rn y \ z 


- 2 

1 + (-Z-) + 
Rn 


- 2 
(— ) 


*T 


1/2 


( 2 ) 


The slope of the shock wave with respect to the free-stream velocity (shock- 
wave angle) is f , and it is determined from the relation 


sin r 



- 2 - 21 - 1/2 
1 + (-Z-) + (— ) 


-n 


^T 


(3) 


In an x = constant plane, the angle a is defined by 

tan a a — — 

B y 

where 


B 


h_ 

Rn. 


(4) 


(5) 


is the ratio of the shock principal radii of curvature at the stagnation 
line. Using eqs. (3) and (4), eq. (2) can also be written as 


e- a - sin T e- + cos T (cos a e- + sin o e-) (6) 

n x y z 

The angles T and a are illustrated in Fig. 2. 

The continuity equation for general three-dimensional flows is auto- 
matically satisfied by introducing a pair of stream functions ¥ and 4> 
such that 


- pV = V4 1 x V4> 


(7) 
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For the special case of axisymmetric flow the second stream function $ is 
identically the circumferential angle $ (where tan <j> = z/y), and hence $ = 
constant planes are | = constant planes. Streamlines are formed by the 
intersection of 4* = constant with <& == constant surfaces. 

For general three-dimensional flows it is assumed here that inside the 
shock layer the surface $ = constant contains generators which are straight 
lines normal to the shock wave, and at the shock wave the velocity vector 
given by the shock relations lies in the 4 = constant surface. Thus the 
$ = constant surface has the correct shape and slope at the shock wave, but 
it is not constrained to satisfy all the flow field equations inside the 
shock layer. For conical flows, this was shown in Reference 5 to be a good 
approximation except near the body surface. 1c is believed that this in- 
accuracy near the body surface has only a small effect on the shape of the 
shock wave itself. It should be noted that for the applications herein it 
is the shape of the shock wave that must be reasonably accurate, not the 
flow field properties inside the shock layer. 

Define a three-dimensional shock-oriented coordinate system £, 3, n 

where £ and B are coordinates along the shock surface with B = $ , , and 

sh 

n is the straight line normal distance inward from the shock wave. With ' 
e- and e^ unit vectors in the directions of £ and B, respectively, 1 
the velocity vector becomes 

V = u e- + v e- + w e- (8) 

s n B 

Since the velocity vector lies in the $ => constant surface and B = 4 1 , then 
w = 0 and 
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( 9 ) 


'-r - A - a 

V = u e- + v e- 
s n 

At the shock wave, e-, e-, and are mutually perpendicular; hence, the 

direction of e which makes w => 0 at the shock is obtained from the shock 
s 

relations as 


e- « cos T e- + sin f ( cos a e- + sin a e-) 
s x y z : 


and thus 


e_ o e- x e- = - sin a e- + cos a e- 
8 s n y z 


( 10 ) 


( 11 ) 


These results may also be obtained from eqs. (37) and (38) in Reference 1 
by setting 6 “ 0 in those equations. The curve n = 0, 8 “ constant on 
the shock wave represents the locus of points whose tangent is the velocity 
component parallel to the shock wave. These curves on the shock surface are 
the same as the simplified streamlines discussed in References: 1 and 2, 
except there the lines were on the body surface instead of the shock wave. 

From eqs. (7) and (10) it follows that 


and 


pu 


1 

h 3n 


( 12 ) 


1 ^ 
pv “ 

h h 9$ 
s 


(13) 


where dS « h d£ is the differential arc length along the 5 coordinate, 
s 

dq « h d|3 is the differential arc length along the 8 coordinate, and 
h g and h are the scale factors or metric coefficients corresponding to 
the £ and 8 coordinates, respectively. 
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Comparing eqs (12) and (13) with those for axisymmetric flows (Ref. 6), 
it is observed that the results here are similar if h replaces the axisym- 
metric body radius and dS = h g d£ is the distance along the axisymmetric 
body. Thus, the assumptions made here, are somewhat like an axisymmetric 
alanogue for the invisicd flow field. 

The value of the first stream function ¥ on the shock wave represents 
the mass flux (per unit g) entering the shock layer. The mass flux per unit 
area is p V sin F and the differential of surface area on the shock is 
hdg h^ d 5, Hence the mass flux per unit $ is 


Y , = p V sin r h h d? (14) 

sh J “ 0= s 

J O ‘ 

For the special case of axisymmetric flow, h = r^, sin r h g d£ = dr g ^ and 
hence the familiar result 


sh 


■ »» V „W 2 


is obtained, The shock standoff distance can be obtained by integrating 
Eq, (12). However, before this integration can be performed, a number of 
subsidiary relations are needed, and these relations are developed below. 
Following the technique used in Reference 1 and 2, equations (10) and 
(11) are used to obtain the transformation operators 


and 


1 3 p 9 , . — / — 3 , . - 9 \ , K . 

7— — - = cos F — — b sm r cos a — — b sm a — (15) 

h g 95 9x \ 9y 9z j 


19 . - 9 ... - 9 

— — = - sin o b cos o — 

h 9g 9y 9z 


( 16 ) 
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The line 3 = constant on the shock wave follows from Eq« (15) as 


3_z 

ll 

n 


1 z 

tan a ® — — 

5 y 


( 17 ) 


which integrates to 


— B — “ 

z a C y 


(18) 


Note that the parameter C is a function of 8 (on the shock wave) . 
Equation (18) can be differentiated to give 

38 C 


and 


3y C'y 


li- 

3z C'z 


(19) 


( 20 ) 


where C* « 


dC 

d8 


Using equs«. (19) and (20) above, substitute 8 into the operator of Eq. 
(16) to obtain 


h ® 


-B+l 

z 


^ [z 2 + B 2 ^ 2B /C 2 ) 1/2 


( 21 ) 


Now along the shock wave (with 8 “ constant); 


dS “ h g d£ » [(dx) 2 + (dy) 2 + (dz) 2 ]- 


1/2 


( 22 ) 
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_ _ _2 —2 - 2 B —2 _ _ o — — 2 ^2 

dS = (dz/z)[z^ + B Z Z za /C^3 [1 + U/\) + Cy/Rj'i') ] 


( 23 ) 


Finally, substitute eqs. (4), (21), and (23) into (14) to obtain 


cr ; B + 1 


sh c 2 (B + 1) 


p v 

r CO 00 


(24) 


In a manner similar to that of References 1 and 2, the parameter C(3) 
is determined by assigning a value of 3 to each line on the shock surface 
at some point hear the stagnation line (see Fig. 3). This is accomplished 
by setting 

£ 2 £ 2 x- (25) 

and 



(2i) 


Thus at x » x- eq. (1) gives 
c 

-2 2 - 
2 6 sin 3 

2— as ■ ■ ■ — » —■ ■■■• ■ 

sin 2 3 + B cos 2 3 

and 

o “2 2 - 

- 2 s cos 3 

7 — B l . n , . 

c 2 - 2 - 

sin 3 + B cos 3 


i 

(27) 


i 


(28) 


Then eq. (24) gives the result (for 3 “ constant) 


sh,e 


as 



(29) 


15 



where 


and 


y _ B (sin 2 3 + B 2 cos 2 B)e 2 y 
S ^ ,e (B 4- i)(sin 2 6 + B gos 2 3) 2 

2 ,(B-l)/2 


C(3) 


B — 
sin B 

cos 3 


2 - - 2 - 
sin 3 + B cos 3 


(30) 


(31) 


With eqs. (3), (18), (27), (28) and (31), sin T near the stagnation line can 
be written as 

f _2 - - 2 2 - -2 2 - - - 2B ^^ 2 

I s [(z/z-) sin 3 + B cos 3 (z/z-) ] 

sin ? » <1 + — 


(sin 2 3 + B cos 2 3)R^ 




(32) 


and using eq. (29) it can also be described as 

2/ (B+l) 


-2 

e 


sin T •=> <1 + 


V 


2B , 


Bi “ 2 S < T * 


+ B 2 COS 2 3 (V , /¥ . -) B+1 
sh sh,s 


(sin 2 3 + B cos 2 3) 


- 1/2 


(33) 


The relation for the stagnation line shock standoff distance A follows 
from eq. (12) as 


A » - 


¥ . 

Sh dV . 1 
o pu 


tl d(V/V sh ) 


P 2 J 


o (uh/f sh ) 


(34) 


Along the stagnation lines both ¥ and V ^ are zero. However, the ratio 
V/'f h approaches a value in the range 0 < 'f/'J'^ f 1 in the limit as the 
stagnation line is approached ^ “ 1 on t ^ ie shock and ¥/¥ 53 0 on 

the body). Also, the ratio uh/’f g ^ is indeterminate along the stagnation 
line, and care must be exercised in evaluating this indeterminancy. The 
velocity component u is obtained from the approximate energy equation used 
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by Maslen (Ref. 4), i.e. 


-2 

u = 


20> 2 - h e ) 


(35) 


where is the static enthalpy aft of normal shock wave and h £ is the 

local static enthalpy. Since h g is a function of the pressure and entropy, 
a Taylor series expansion about the position a streamline crosses the shock 
wave yields 


h e (€,M) “ h sh <M) + 


/ ah 

e 

[3P' 


sh 


tpcc.e,'?) - p 8h (M)] 


(36) 


For isentropic flow along a streamline, the Bernoulli and energy equations 
give 


f Sh 
€ 

i3p 


'sh p sh (e,\i>) 


(37) 


With p = p(p,h e ) a Taylor series expansion about the normal shock position 
produces 




JL_ 

P, 


+ t 


3(l/p) - 

8p 




P 2 ] 


+ [ 


5 (1/p) - 


9h 


th sh (3,^) 


h e> 


(38) 


In order to determine h g from eq . (36), an expression for the pressure 
across the shock layers is needed. The momentum equation normal to a stream- 
line is 


i£ 

9h 



(39) 


As in Reference 4, the streamlines are assumed to run parallel to the shock 
wave in a 6 = constant plane; therefore 
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V - u 


and 


R = R 


Then substituting these results and eq . (12) into eq. (39), the following 
results is obtained 



(40) 


Finally assume u/(hR) is constant across the shock layer and equal to its 
value on the shock wave. Then eq . (40) becomes* 


ip. 

34 * 


u sh 

hR 


(41) 


Equation (41) integrates to 


where 


P(C,B,Y) = P sh U,0) + 


Rh sh 



(42) 


(43) 


is the shock radius of curvature in £ direction. The distinction between 

P .($,’*') and p ,(£,B) should be noted in the sketch below, 
sn sn 


As noted in Ref. 4, the assumptions leading to eq . (41) introduce compen- 
sating errors, and the resulting pressure equation is quite accurate 
except near the stagnation line. 
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shock 



Next, define 

sin 2 r = 1 - G 2 (44) 

p sh 5 P 2 " Ap ° 2 (45) 

h sh 5 h 2 “ Ah ° 2 (46) 

where near the stagnation line 

G 2 << 1, Ap G 2 << p£> Ah G 2 << h^ 

and Ap and Ah depend only on the fluid properties aft of the normal shock 
wave and are obtained from the shock-wave relations for a perfect gas and 
equilibrium air in Appendix A. Now, substitute eqs. (36), (37), (38), (42), 
(45), and (46) into (35), expanding and neglecting the higher order terms; 
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the following result is obtained 


u 2 (e J.v) 


2Ap 


2 - - Ahp„ 2 - U sh ^sh V 

G‘ U,S) + - l) g 2 (B,v) 4- (1 - J-) 


Ap Rh 


sh J 


(47) 


and hence 


h 2 u 2 (g,B,¥) = 2Ap. 


r r2„2 - 




sh 


Ahp 


r2„2 r.v u 


h G (c,B) , , M 2 h G (B,y) , sh h 1 , V . 

n • ^ 'A- . o * r _ A ~ ' ■*- “ \U / 


r 


Ap 


T 


R ¥ . A P 
sh 




sh 


sh sh 

Expanding eq. (32), neglecting the; higher order terms, and comparing 
the result with eq. (44); the following relation is obtained 


(48) 


g 2 U,6) 


-2 

6 


- - 2 2 - -2 2 - - - 1 2B 
(z/z-) sin 8 + B cos 8 (z/z-) 
c e J 


R^ (sin 2 8 + B cos 2 8) 


(49) 


In a similar fashion eq. (33) yields 

2/(B+l) 


G 2 (M sh ) 


-2 

e 


sin 2 I (y^;) 


+ * 2 cos2 * ( 


2B/(B+1) 


R^ (sin 2 8 + B cos 2 8) 


(50) 


This equation is important because it establishes the functional form for 


G 2 (B,'y). Replacing Y , with W, it gives 


G 2 (8 »H') 


-2 

e 


sh 

2/fB+l) 2B/(B+1) 

sin 8 ( V F/'R . -) + B cos 8 (¥/¥ . -) 

. sn > € sh | £ 


R 2 (sin 2 8 + B cos 2 8) 


(51) 


Using eqs. (21), (23), (24), (32), (49), and (50) the following results 
(which are required to integrate eq. (34)) are obtained in the limit as the 
stagnation line is approched. 
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£im h 2 G 2 g»e) 

* 2 , 


(B+l) _ 
2 ,. 2-2 


( 52 ) 


For every position away from the stagnation line, the ratio is unity 

on the shock wave and aero on the body. Thus, in the limit as the stagnation 
line is approached (z -*■ 0) , 


0 < H'/'i' < 1 

sh - 


although both y 0 and ¥ ^ -> 0. The remaining terms approach the follow- 
ing limits: 

Urn k 2 G 2 (&,^) r&im h 2 G 2 (LI) ; riim G 2 (g,V) . 

M L ' M 4 *° <= 2 <^ 


(B+l) Aim j- 

pVi^ w 


sin 2 8 (^/^ sh “) 


2yf+l) _ 2 2 _ 


+ B cos 8 (¥/¥: , -) 
sn » s 


2B/(B+1) 


- 2AB+1) 2 7 2B/(Bd-l) 

si “ s <f S h rt sh,; ) + B cos » (< ' 8 h /T S h,; ) 




(53) 


where 


A = < 


2B 

B+l 

2 

B+l 


for 0 < B < 1 


for 1 < B < » 


(54) 


Since u *= cos T (from the shock wave relations) , 
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Jlim U sh ^ 

if. 

sh 


l im 3 (sin F) 


iim 


z+0 V 


sh 3s 


h V 


z-*0 f 


sh 


3 (sin F ) \ 32 


3z 


3S 


Aim (B+l) 
z ->0 5 2 

P “ 


. *3 -2B 

-2 , B z 
z + 


li 


-2. , B 2 5 2B1 

Z + 


c 2 


: -2 -2 -2Bi2 

i + ^ + B 2 


-2 -2 

*T C 


( 55 ) 


- 


\ 


(B+1)B 

P. 


(B+l) 

P- ^ 


for 


for 


0 < B « 1 


1 < B < “ 


(56) 


Finally, substituting eqs . (48) and (51)-(56) into eq. (34) the follow- 
ing equation is obtained for the shock standoff distance 


where 


and, 


®T 1 


p V p -11/2 . 

( (— ) ■ - 1 -- 

2Ap VJ (B+l) * 


/ 2 
p ® V ® B 


D = < 


Ap B+l 
V 2 . 

00 00 J. 

Ap B+l 


dn 


o [1 + D(1 -tO + Eq A ] 1/2 


' sh 


for 


for 


0 < B < 1 


1 < B < «° 


g a — — P ”--” 

V 2 p » Ap 
00 


- 1 


(57) 


(58) 


(59) 


(60) 


22 



Expression for ^ 2 — - and are given in Appendix A for both an ideal 

p V V 

CO oo oo 

gas and equilibrium air. For an ideal gas the results from Appendix A are 


and 



2 

Y+l 


E = 


0 & - 1> 2 

2M 2 (1 + m 2 ) 

« 2 00 


In general, eq. (57) cannot be integrated in closed form. However, for 
the special case of B <= 1 (spherical shock at stagnation line) it can be 
integrated to yield the same result given by Maslen in Ref. 4, i.e. 


A a 

r p V 2 p - 

oo oo oo. 

1/2 

1 



R 

L 2Ap p 2 j 


1/2 

p v 2 

OO CO 

1/2 - 




(1 + E) ' + 

I 1+ 2 p 



(61) 


which compares well with more exact numerical results and experimental data. 


Shock Shape Downstream of Stagnation Region 

In reference 1 a method is developed for calculating the inviscid sur- 
face streamlines from a prescribed pressure distribution assuming "normal 
shock" entropy at the edge of the boundary layer. Before the entropy at the 
edge of the boundary layer can be calculated, the shock-wave shape associated 
with each inviscid surface streamline must be determined. Normally, the : 
entire shock shape in the subsonic-transonic region would have to be deter- 
mined iteratively. However, it is assumed here that for a given surface 
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pressure distribution the shock shape associated with each inviscid surface 
streamline can be calculated independently of the others. This' procedure 
highly simplifies the calculations and is consistent with the accuracy of 
the assumptions used in Reference 1. 

For the region downstream of the stagnation line, it is assumed that 
3 = 3, i.e. a 3 = constant plane is assumed to intersect the body surface 

on a line of 3 “ constant which is an inviscid surface streamline. With 
this assumption, eq . (42) can be applied to the surface to give 

u . (?,e) v . 

p h U|i*) a p ah (c,I) r (63) 

b sn R h 

Then using eq . (43) 


ii V V V - vv - 

sh sh = _ °° sh 3 (sin T) _ <*» sh 5 (sin T) 

- - - h 3S 

R h h 3S 


(64) 


where h is the scale factor for the coordinate 3 on the body surface and 
S is distance along an inviscid surface streamline. 

Now substitute eq . (64) into (63) to obtain 


9(sin f) 
3S 


(p b " P sh^ V V 


“ sh 


(65) 


This equation can be integrated numerically along an inviscid surface stream- 
line to determine the value of sin F (where T is the shock angle) 
associated with each position. The value of for each position can be 

calculated numerically from eq. (14) in the following form 

hdS (66) 

O O 


f , = p v 

sh » » 


sin f hdS - p V 

00 OS 


sin f 
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Thus for each position along an inviscid surface streamline the correspond- 
ing values of r and Y ^ are known. As will be shown below, the value of 
? (shock-wave angle) corresponding to a given value of ¥ ^ is all that is 
needed to determine the entropy at the edge of the boundary layer. The sur- 
face pressure distribution used in eq, (65) must be modified so that 
Pb = p^ at the stagnation line instead of the stagnation pressure. For a 
given surface pressure (p ) distribution which yields p = p ■ at the 

£ £ Stclg 

stagnation line, this can be accomplished by simply using 

(P 2 “ O 

p, = P + (p “ P ) 7 7 (67) 

r b w *e <*>' (p - p ) 

*s °° 

where p is the pressure to be used in eq. (65) only. 

Entropy at the Edge of the Boundary Layer 

The mass flux within the boundary layer along an inviscid surface . 
streamline is 

p u h df$ ( 6 - 6 *) 
e e 

where p g u e (6 - 6 *) is the mass flux per unit width and hdS is the 
spacing between inviscid surface streamlines. 
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Therefore the mass flux per unit 0 is 




B.L. 


p u h (6 - 6*) 
e e 


( 68 ) 


where 6 is the boundary layer thickness and 6* is the displacement 
thickness. Now the streamline at the edge of the boundary layer passed 
through the bow shock wave at the position where 


'v = y 

sh B.L. 


(69) 


or, in other words, where the mass flux (per unit 8) inside the boundary 
layer is equal to the mass flux (per unit 8) entering the shock wave 


(see sketch below). 



Then for the value of ¥ , given by eq. (69), the shock angle T is 

sn 

determined by the method developed in the previous section. After f is 


26 



determined, it is a relatively simple matter to calculate the entropy for 
that streamline from the shock relations. 

The procedure described above must generally be performed iteratively 
because the properties on the right side of eq. (68) require both the pres- 
sure and entropy at the edge of the boundary layer. Since the entropy is 
not known initially, a value of normal shock entropy (or value for last 
station along that surface streamline) could be used to start the iteration. 
Convergence is quite fast and generally requires less than four iterations. 


Calculation of 6 - 6* 


Laminar Boundary Layer 


In reference 1 a method was given for calculating the momentum thick- 
ness 0^ based on applying the axisymmetric analogue to Beckwith and 
Cohen's boundary layer results (ref. 8). Following the technique of 
reference 8, it follows that 


6 - <5* 


. m 


e > : s 


(70) 


and from Table 1 in reference 8 

8 = .4696 and 6 = 1.2168 

tr tr 

In addition, it is assumed that <5*^ » 5.2//J, corresponding to the factor 
5.2 used to define the edge of the boundary layer in the Blasius solution 
for a flat plate. Using these results, eq . (70) becomes 
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(71) 


6 - 6 * - 5,27 0 

m 


for the laminar boundary layer. 
Turbulent Boundary Layer 


In reference 1 the turbulent boundary layer theory of Reshotko and 
Tucker (ref. 9) was used. Following the method of Ref. 9, 


5 - < 5 * - 6 (~ 

m t) 


- H) 


m 


Vef - « 


(72) 


Then with a 1/7 power law for the transformed incompressible velocity 
profile, it follows that 

. (7 + I )(-6 + 2) , . 10 3 
tr 


Using this result in eq, (72) the following relation for 6 - 6 * is ob- 
tained 


6 - 6 * = 6 (10.3' - H) (73) 

m 

for the turbulent boundary layer. The form factor H.= 6*/0 m is calculated 
by the method described in reference 1 . 

Transition Region 

For the transition region, the weighting function w^ used in reference 
1 is used to give an average value of 6 - 6 * as 
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(74) 


6 - 4* - (1 - w f >5.27 + w f 6^, (10,3 - H) 


where w^ = 0 for laminar and w^ = 1 for turbulent flows. 


Method for Calculation 


The technique described in this report is coupled with that of reference 
1 to calculate the heat transfer rate with variable entropy at the edge of 
the boundary layer. Appendix A describes the method used to determine B 
and R^, on the shock wave at the stagnation line which correspond to the 
body shape (B and R^ for the body differ from the corresponding quantities 
on the shock wave) . The changes needed for the computer program listed in 
reference 3 to be modified to include variable entropy at the edge of the 
boundary layer are given in Appendix C. 
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RESULTS AND DISCUSSION 


Table 1 shows a comparison of the stagnation line standoff distances 
computed by the present method, eq. (57), and Maslen's method (ref. 5). 
Maslen obtained different results depending on the direction that the stag- 
nation line was approached.. Thus the r.m.s. value is given with the r.m.s. 
deviation from that value. This table shows that the present method com- 
pares well with Maslen's for the stagnation standoff distance. 

Figure 4 gives the variation of A/R^ with B for M m =» 2, 5, and 10. 

As expected, the differences between the M^ = 5 and 10 results are small. 

Blunted 15° Half-Angle Circular Cone 

Laminar heating rates were calculated along the windward plane of 

symmetry of a spherically blunted (R^ =■ 0.375 in.) 15° half-angle circular 

cone at M » 10.6, a = 20°, and Re . = 3.75 x lo\ Since the nose radius 
OO »,N 

is a significant parameter for variable entropy effects, the Reynolds number 
Re^ N is based on the nose radius of the body. In order to compare the 
numerical results with Cleary's tabulated experimental data (Ref. 10), the 

free-stream properties corresponding to wind-tunnel stagnation properties 

5 2 * 

of p = 1.728 x 10 lb/ft and T *» 2000°R were calculated to be 
s s 

P m - 2.6614 lb/ft 2 , T w = 89. 971°R, and - 4928.1 ft/sec (M w - 10.6). 

Also, it was determined that a value of 4 = h /H « 0.251 corresponds to 

T w /T g = 0.270. Gas imperfections in the wind-tunnel stagnation properties 
were taken into account in determining the free-stream conditions; however, 
the laminar heating rates were calculated using a perfect gas with y « 1.4. 
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The experimental heating rates in Ref, 10 are ratioed to the calculated 

2 

value of q = 35.94 Btu/ft -sec in this report; 
n w , s r 

Figures 5 and .6 compare the calculated ratio of local to stagnation- 
point heat-transfer rate along the windward plane of symmetry with the ex- 
perimental data. In each figure numerical results are presented for both 
variable entropy and normal-shock entropy along the edge of the boundary 
layer. Figure 5' gives the. heating-rate ratios using inviscid surface 
streamlines calculated from modified Newtonian pressures. When the vari- 
able entropy effect is included, these streamlines can differ significantly 
from those calcual ted using normal-shock entropy. Thus, there are essen- 
tially two factors involved in the variable entropy effects on the heating 
rate - one- is the scale factor associated with the streamlines and the 
second is the flow-field properties at the edge of the boundary layer. The 
first factor tends to decrease the heating rate when the scale factor de- 
creases, whereas the second factor tends to increase the heating rate where 
variable entropy effects are significant. In Figure 5 the heating rates 
calculated with variable entropy effects are slightly below those calcualted 
with normal-shock entropy. Upon examining the numerical results it was 
found that the variable entropy effect reduced the streamline scale factor 
h along the windward plane of symmetry which, in turn, reduced the heating- 
rates slightly more than the increase in heating-rates due to the variable 
entropy effect on the flow-field properties at the edge of the boundary 
layer. 

Figure- 6 shows' the calculated heating*-rate ratios- using the simplified 
streamline patterns discussed in Ref. 1. The simplified streamlines depend 
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on the body geometry and angle' of attack only; hence-,' variable' entropy has 
no effect on the- geometry of the simplified streamlines and their .corre- 
sponding scale factors. As a result, the variable' entropy effect increased 
the heating-rate ratios in Figure 1 . 6 about 15% higher ' than those calculated 
using normal-shock entropy-. The experimental data fell between the two 
theoretical solutions » 

The heating rates for the cases presented in' Figures 5 and 6 were 
also calculated using tangent wedge and tangent' cone pressures along the 
conical afterbody in lieu of that given by eqv~( 67) . ' However , no discern- 
ible differences could be detected when the results' were plotted on 
Figures 5 and 6', and therefore they are not shown. 

Blunted 2:1 Elliptical Cone 

""Laminar heating rates were also calculated on a blunted 2:1 elliptic 

cal cone at . a = 15°, 30°, and 60°; M = 10; and Re ® -2.61 x lo\ 

99 *00 * °°|N 

4 5 

8,39 x 10 , and 2,02 x 10 . Figure’ 7 illustrates the' geometry of the 
model used for the experimental'heat-transfer- data- presented in Ref. 11. 

The pertinent flow-field properties required as input data' for each case 
calculated in this report are presented in Table 2. All lengths and 
Reynolds numbers. (Re^ N ) are referenced to the nose radius in the plane of 
the major axis (Rj^ = 1.0 in.). As noted in Ref. 2, modified Newtonian 
pressures are not accurate enough to- give the proper spreading of the 
streamlines (and hence the correct scale factor) for this body shape. 
Therefore, the simplified streamlines discussed in Ref. 1, which depend on 


32 



body geometry and angle of attack only, are used in. each case presented 
herein. 

Figure 8 shows the heating-rate ratio along the windward plane of 
symmetry for Re^ N =» 8.39 x 10^ and a ® 15°, 30% and 60% Downstream of 
the nose the variable entropy effect increased the heating rate ratio about 
10% for a 15°, 20% for a « 30% and 30% for a « 60’’’ over that calculated 
using normal-shock entropy. These percentages tended to remain constant in 
the downstream region. Heating rates calculated with variable entropy 
agree much better with the experimental data than those calculated with 
normal-shock entropy. The effects' of variable entropy on the heating rates 
begin closer to the nose for the larger angles of attack cases. 

The effect of Reynolds number on the laminar heating-rate ratio with 
variable entropy is depicted in Figure 9 for the windward plane of symmetry 
on the blunted 2:1 elliptical cone at a « 30°. Since the boundary layer 
thickness increases as the Reynolds' number decreases, the effect of vari- 
able entropy on the heating rate begins closer to the nose for the smaller 
Reynolds numbers. This effect can also be observed in Figure 10 which 
shows the variation of the shock angle corresponding' to the position where 
the streamline at the' edge of the boundary layer crossed the shock wave. 

This shock angle decreases downstream of the nose with the smallest Reynolds 
number case decreasing the fastest. For x/R^ 2 10 the shock-wave angle 
for all three Reynolds number cases approaches the' straight portion of the 
shock wave surrounding the conical afterbody in the windward plane of 
symmetry. 
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Figures 11 and 12 show the circumferential" variation' of the . heat „ 
transfer coefficient’ on the blunted 2:1- elliptical: cone at a “ 30® and 
60°, respectively, for the axial positions- of 4.7’ and 9.7. The 

heat transfer- coefficients calculated" with- variable’ entropy effects agree 
much closer with’ the" experimental -data -than" those’ with- normal-shock entropy 
on . the -windward surface. However; very"little difference- was found for the 
leeward region. Figure” 13’ shows ' the 'drcumferential variation in the 
shock-wave angle” corresponding to the streamline at the' edge of the bound- 
ary layer at’ axial stations of x/R^ «"2.2, 4.7, and 9,7, These results 
show that the shock angle decreases going” away from the windward plane 
(<t>“0) , reaches a minimum value near " $• “ 100° , and then increases as the 
leeward plane (<£=180°) is approached. 

Space Shuttle Orbiter 

Figure 14 shows a typical’ space shuttle" orbiter" (without a canopy) 
on , which computations have been performed,' The analytical description is 
similar to that given in Ref v 12; it employs- eTlipticaT’ cross-sectional 
shapes with different ellipticity- ratios - on windward and” leeward sides. 
Cubic polynomlnals are used to define 'the plan' and thickness distributions. 
Based on a length of L- = 1ft. , the region for' x/L <0,05 is a spheri- 
cally blunted 15® cone whose nose radius is R^/L *»' 0,01512, For 
0,05 < x/L < 1.0, the body- shape is the- same” as' that' in Ref. 12 where the 
body radius f(x,<j>) is given by 

- 1/2 

f (x, «J>) = yfcos 2 ■•+ (y/z) 2 sin 2 <j>] (75) 
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where 


) = a + a.x* + a 0 x* 2 + a.-X*'* (76) 

( o 1 2 3 

z J 

x* * x - x^ (77) 

and the coefficients for the equations above are given in Table 3 for five 
body segments. This description of the body produces continuous body 
slopes but discontinuous radii of curvature. 

As noted in Ref. 2, since the windward side of this vehicle is nearly 
flat the modified Newtonian pressure distribution does not give an accurate 
description of the inviscid surface streamlines. Therefore, the results 
presented here calculated the streamlines by the simplified method given 
in Ref. 1, which depends only on the body geometry and angle of attack. 

The following properties' were used for all the calculations: 

p w - 2. 2766 lb/ft 2 , T^ -■ 87,31°R, « 4639.42 ft/sec, (M^ - 10.13, 

Re = 1.5 x 10 4 ) , X> * 0.2501, and v “ 1»4, 

Figure 15 shows the laminar' heating-rate ratio along the windward 
plane for a = 15°, 30°, and 60®. The trend here is somewhat similar to 
the blunted 2:1 elliptical cone in that the variable entropy effects in- 
crease the heating-rate ratio more at the larger angles of attack. Figures 
16, 17, and 18 illustrate the circumferential heating-rate ratio at. 
x/L =>• 0,1, 0,3', and 0.5, respectively. These figures show that the vari- 
able entropy effects on the heating are significant on the windward side 

but insignificant on the leeward side. Figures 17 and 18 also show that 
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the, heating-rate ratio reaches a maximum. off the windward plane of symmetry 
for x./L = .0.3 and 0.5. 

Heating-rate ratios for laminar, transitional, and turbulent, heating 
along the windward plane of symmetry are presented in Figures 19 (a) & (b) for 
a = 15° and 30° , respectively » The transition region was arbitrarily chosen, 
to be. 0,55 < x/L ,< 0,80, These figures show that variable entropy effects 
on the heating are much more pronounced in the turbulent region than the, 
laminar religion. For a =30°, the increase in the calculated heating-rate 
ratio with variable entropy over that for normal-shock entropy is about 10% 
for the, laminar region whereas it is about 50% in the turbulent region. 


CONCLUDING REMARKS 

A relatively simple method is presented for modifying the heat trans- 
fer method developed in, References 1-3 to include the effects of variable 
entropy at the edge of the boundary layer » An approximate shock-wave 
shape corresponding to. each inviscid surface streamline is calculated 
using a modified form of Maslen's method for inviscid axisymmetric flows 
(Ref, 4), The position where the streamline at the edge of the boundary 
layer crosses the shock wave is determined by equating the mass flux pass- 
ing through the shock wave to , that , inside the boundary layer. The slope 
of the shock wave at that position determines the entropy of. that stream- 
line, The approximations used in this technique allow the shock-wave shape 
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and surface heating rates corresponding to each inviseid surface streamline 
to be calculated independent of the other streamlines,, 

The approximate inviseid solution used to generate the shock-wave 
shape is based on.Maslen’s observation that most of the mass flow inside 
the shock layer is concentrated near the shock,, Hence, it is assumed that 
across the shock layer the streamlines are essentially parallel to the 
shock and have the same direction as the velocity vector just aft of the 
shock wave. This assumption is somewhat like an axisymmetric analog for 
the inviseid flow field,, Although the approximations involved are inac-^ 
curate near the surface, they appear to have little effect on the shape of 
the shock wave itself. It is the shape of the shock wave that is needed 
here and not a detailed structure of the flow-field properties across the 
shock layer. 

The shock standoff distances given by the present method compare well 
with those computed by Maslen's more elaborate asymmetric method (Ref. 5). 
In order to make the computation of the shock shape downstream of the nose, 

a direct one, an approximate surface pressure distribution (eq„ (67)) was 

used. This distribution was found to yield accurate shock shapes for a 
sphere, but it has not been checked for other body shapes. If the surface 
pressure distribution was not specified, the shock shape solution would 
have to be solved by an iterative scheme. 

The calculated laminar heating rates on blunted circular and ellip- 
tical cones showed that the effects of variable entropy were significant 
on the windward side but insignificant on the leeward side. At both 

a “ 30° and 60° the calculated heating rates on the blunted 2:1 elliptical 
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cone were very close to the experimental data with variable entropy, where- 
as the heating rates with normal-shock entropy were significantly lower 
than the experimental results on the windward side. It was found that 
variable entropy effects started nearer the nose for those cases calculated 
at the lower Reynolds number. The percentage increase in laminar heating 
rates due to variable entropy effects was also higher at the larger angles 
of attack. 

The calculated heating rates on a , typical space shuttle orbiter 
showed . that variable entropy effects were significant over a large portion 
of windward plane of symmetry. In addition, variable entropy effects in- 
creased the calculated turbulent heating rates much more than it did for 
the laminar heating rates. However, comparisons with experimental tur- 
bulent heating rates are. needed to assess the accuracy of the present 
method. 

The method presented here requires comparatively small computational 
time. A typical case requires only a few seconds per streamline on the 
IBM 370/165 computer. 
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APPENDIX A 


DERIVATION OF Ap AND Ah 


The oblique shock-wave relations for an equilibrium or ideal gas are: 


p V sin r = p , V , 
» 50 sh sh 


(Al) 


2 2 - 

p + p V sin T = p , + p V sin TV, 

r CO OO 00 sh 00 00 sh 


(A2) 


h + 


„2 . 2 - 2 

V sin T V 

” v , sh 

2 h sh + 2 


(A3) 


2-2 2 2 
Substitute sin T = 1 - G , p gh = P 2 “ ApG , h gh = h.^ - Ah G (where 

2 

G << 1), and eq. (38) into the above equations. Neglecting terms of -order 


higher than G , the following results are obtained: 

(p./p, - (1 - P>> - 1]. 

^£_ = _ 


and 


P V 




3 v 2 


<30 00 

4h % ' “ [ ~ 9 ? 




2 v p 


3p 


i* pf t ^ + * t ^ 2 


(A4) 


(A5) 


For the special case of an ideal gas. 


i . (izi) ^ 

P Y P 


(A6) 


aq/p) 

8p 


1_ 

PP 


(A7) 
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and 


iUM = (ini) 

ah e tp ; 


(A8) 


Then using the normal shock relations to obtain expressions for &2* 

and h^ and substituting eqs. (A6)-(A8) into eq. (A4) yields 


Ae. 



CO 


2 

Y+l 


(A?) 


Further, combining eqs. (A5) and (A9) and using eqs. (A6)-(A8) yields the 
following result for the ideal gas 


Ah ^2 P °° V « . 
= „2 P Ap 


or - i). 


2M^ (1 + ^ M 2 ) 


(A10) 


For equilibrium air, the correlation formulae of Reference 7 are also 
used here to be consistent with the approach taken in Ref. 1. The expression 
for the reciprocal of the density is (Ref. 7) 


x 1 - 1.0477[1 - (h e /h E )- 6123 ] 


7.344 x 10 5 (yfjy) 


.965 


(All) 


where 


h_ = 2.119 x 10 3 ft 2 /sec 2 


(A12) 


From the above equation the following expressions are obtained 


(1/p) -I .965 

‘»P 2 ' P2»2 


(A13) 
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and 


. 6123 

.8(1/0), U.0477)(.6123)(h 2 /h E ) 

3h e 2 " -5 ”2 ' 965 


h 2 7.344 x 10 


( — — ) 
''2117“' 


(14) 


The normal shock poperties p 2> p 2 and h 2 are calculated by the iterative 
technique described in Reference 1. Note that any other equilibrium gas 
model could also be used in lieu of the one described above. 
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APPENDIX B 


RELATIONSHIP BETWEEN SHOCK AND BODY RADII 
OF CURVATURE AT STAGNATION LINE 

The expression for A/R^ given by eq. (57) requires the shock radii of 
curvature (R^, and either B or Rjj) to be known at the stagnation line before 
the standoff distance A can be calculated. For blunt-nosed bodies, this 
would generally require an iterative solution for the entire subsonic- 
transonic region before R^, and B could be determined. On the other 
hand, these quantities can be calculated if a surface pressure distribution, 
consistent with eq. (42), is known in the vicinity of the stagnation line. 
For the analysis herein, a surface pressure distribution similar to modified 
Newtonian pressures is used for calculating the shock wave shape. This 
pressure distribution is 


p b -■ (p 2 - p«,> sin 2 r b + Po 


(Bl) 


where is the inclination of the body surface with respect to the free- 

stream velocity (r^ « tt/2 at stagnation point). Note that this equation is 
a modification of modified Newtonian pressures in order to give p^ « p£ at 
the stagnation point, which is required to be consistent with eq. (42). 
Equation (42) also gives the surface pressure (by substituting V = 0) as 


p b ’ + 


T sh u sh (;,B) 
R h 


(B2) 


Equating the right sides of eqs. (Bl) and (B2) , there results 
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(B3) 



Equation (B3) will give the desired relationship between shock and body, 
radii of curvature] however, the body points (y^» z^) corresponding to 
shock points (y, z) must be established first. 

Let be the position vector for a shock point and the corre- 
sponding body point in the wind-oriented coordinate system, as shown in the 



Then, with e the unit (outer) normal to the shock wave and A the stand- 
n 

off distance, the body point determined by a straight line normal from the 


shock wave is 
Let 





(B4) 

(B5) 
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R «* x e_ + y e~ + z e- 
s x J y z 

(B6) 

Now substitute eqs. 

(B5) , (B6) , and (3) into (B4) . The result gives 



y^ = y - A cos f cos c 

(B7) 


z, <=> z - A cos F sin a 
o 

(B8) 

Near the stagnation 

line, eqs. (4) and (5) give 



cos ? cos a “ y/Sn 

(B9) 

and 

cos f sin o “ z/R^, 

(BIO) 

Substitute eqs. (B9) 

and (BIO) into (B7) and (B8) to obtain 



" * y ( 1 - ~) 

R 11 

(Bll) 


z “ z ( 1 - ~) 

(B12) 


«1 

Now near the stagnation line the body is represented by a portion of an 
ellipsoid, in a manner similar to the shock wave. Therefore, following 
the form of eq. (4), 


sin 2 ? <= 

b 


1 + (■ 


*b 2 


Z K 2 

) + (— ) 


R 1 1 


,b 


^T,b 


-1 


(B13) 


and substituting eqs . (Bll) and (B12) , 
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I - 2 2 -2 2\-l 

sin 2 ? b - 1 + r ^ L ~ (1 - A/R n ) + rf— (1 " A/5,.) (B14) 


•Next, substitute eqs. (B14), (45), (49), (21), (24), (43) into eq. (B3) to 


obtain 


r. z 2 a - A/A ) 2 y 2 (l - A/R U ) 2 .T_ - 
(P 2 “ P J 1 + 2 + 2 1 + p » 


p 2 - Ap ^2 + -a 


-2 -3-21 

i_ + Bjl 

r\ 1 rv I 


, 12 P. V » 

^ J (B+l) fl + L + IS! 


(B15) 


Expand this equation and neglect all terms of magnitude higher than z 
and/or y 2 . Then the resultant must be independent of z and y; hence, 
the coefficients of z and y must each be zero. Equating the coefficients 


of z in the resulting equation gives 


P 2 P °° 1 _ A. 

,p, ^JL 5 

a p + JL 

p V 2 B+l 


(B16) 


Similarly, equating the coefficients of y gives 

\ _ Ml 2 A£_ + -i_ 

-2 L rJ p V 2 B+l 

B T L » e® 

t £3 '■ ' 

B 2 x _ A_ 2 Ap + _J_ 
^ 0. V* B+l 


(B17) 
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where 


B = 


= iiL 


R 


n, b 


(B18) 


In general, eqs . (B17) and (57) must be solved iteratively to find A/R^ 
and B for a given, value of B and Ap/ (p^V^) „ Then eq„ (B16) can be 
used to calculate R^, for a given R^ ^ Note that for an ideal gas, eq. 
(A9) gives. 


Ar 


p 



00 00 


2 

y + 1 ‘ 
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APPENDIX C 


DESCRIPTION OF COMPUTER PROGRAM CHANGES 

This appendix describes the changes necessary to modify the computer 
program listed in Ref. 3 to include the effects of variable entropy at the 
edge of the boundary layer. Changes are required for the MAIN program and 
function RUNGE subprogram in Ref. 3, and two new subroutines SHPROP and! 
SHOFF must be added to the computer program listed in Ref. 3. A description 
of the additional input and output parameters is also given below. 

Main Program 

The main program reads part of the input data, calculates the initial 
data for each streamline on the ^-circle around the stagnation point, and 
then calculates and prints the heating rates and other pertinent data along 
each streamline (independently of the other streamlines). The program 
variables involved with the changes to the main program are listed below, 
followed by a listing of the main program. 

B. ratio of body principal radii of curvature at stagnation point. 

BB ratio of shock principal radii- of curvature at stagnation line;, B 

DOR ratio of shock standoff distance to shock radius of curvature,. 

A/^ 

F(8) sin f (£,&) for RP = 0 

F(9) V . for RP = 0 
sh 

F (10) sin F (f,B) for RP > 0 
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F(ll) a , for KP > 0 ? 

sh 

GBK local. shock angle, f (£,f) 

GBL shock angle where streamline at boundary layer edge crossed the 

" a '” c; ‘ vr "s hock , F (&,¥) 

G(r,L)" ' 4 , g ^ at position L along streamline (slug/sec) 

G(2‘,L) ~ shock angle F (5,3) at position L along streamline 
G(3,L) ' Dr/DV (sec/siug) 


HSHK 

KE 


NE 

PB 

PSH 

PSHK 

PSS 

RHOSHK 

RPER 

RTB 

SPL 

Y(I) 

YB 

2B 


2 2 

enthalpy aft of local shock angle, h g ^ (ft /sec ) 

input parameter; KE => 0 for normal-shock entropy, KE = 2 for 

variable entropy 

number; of differential equations to integrate minus KE 

2 

surface pressure for shock solution only, P^ (lb/ft ) 

pressure aft of shock where streamline at boundary layer edge 

2 

crossed shock, p g ^(lb/ft ) 

2 

pressure aft of shock at edge of local shock layer, p g ^ (lb/ft ) 
stream function, f (slug/sec) 

density aft of shock at edge of local shock layer, p g ^(slug/ft ) 
body radius of curvature, R^, b (ft) 
shock radius of curvature, (ft) 

sine of shock angle where streamline at boundary layer edge 

crossed shock, sin T (g , 'K ) 

f s 

F(I) DS, where I = 1, 2, „ „ „ , NR 
o 

y 

z 
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The listing of the main program and the changes in function RUNGE are 
listed on the following pages. A line drawn through a statement indicates 
that statement in the original program should be deleted. A rectangle 
drawn around statements indicates that these are new statements to be added 
to the program. 
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ISTG IS AN INDICATOR VARIABLE 

IF ISTG=-1 => NEGATIVE ANGLE OF ATTACK 
IF ISTG=G => STAGNATION POINT NOT FOUND 
IF I STG-1 => STAGNATION POINT FOUND 
CALL STGPR (Pi N* TIN f VIN*RHOlN» EMIN, PS*RHOS*HS*KG) 
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S PH = S INI Yl 2 ) J 
C PH = COS{ Y (2) i 
SGM » SINIGKJ 
CGM = CQStGM) 
STH * SIKIYI5I) 



CTH = COS ? Y ( 5 )} 

SDLP SINCDLPHJ 
COLP = COS ( DLPH ) 

FiU = CTH*CGM 

F 12 ) - { STH*CDLP~CTH*SDLP*SGM5 /FF 
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REMAINING PART OF RUNGE; IS UNCHANGED 





Subroutine SHPROP 


This subroutine calculates the fluid properties aft of an oblique 
shock wave for a perfect gas or equilibrium air„ A call to subroutine 
SHPROP has the form 

CALL SHPROP (SG, PSH, RHOSH, HSH, KG) 

where the arguments are 

SG sine of the shock angle (sin T) 

2 

PSH pressure aft of shock, p ^ (lb/ft ) 

3 

RHOSH density aft of shock, p^. (alug/ft ) 

2 2 

HSH enthalpy aft of shock, h ^ (ft /sec ) 

KG indicator variable; KG = .0 for perfect gas, KG ■ 1 for equilibrium 

air 

Other program variables are defined in Refo 3. A listing of subroutine 
SHPROP appears on the following page„ 
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END 



Subroutine SHOFF 


This subroutine calculates the shock standoff distance at the stag- 
nation line and also the shock principal radii of curvature from the body 
radii of curvature at the stagnation point, A call. to subroutine SHOFF 1 
has the form 

CALL SHOFF ( B, RPER, P2, RH02, H2, BB, RTB, DOR, KG) 
where the arguments are 


B 

ratio of body principal radii of curvature 


RPER 

Rrj, , (body radius of curvature) 


P2 

2 

pressure aft of normal shock, p 2 (lb/ft ) 


RH02 

3 

density aft of normal shock, p 2 (slug/ft ) 


H2 

2 2 

enthalpy aft of normal shock, h 2 (ft /sec ) 


BB 

ratio of shock principal radii of curvature, B 


RTB 

shock radius of curvature, R^ 


DOR 

ratio of shock standoff distance to shock radius 

of curvature, A/Rj, 

KG 

indicator variable; KG =0 for perfect gas, KG = 

1 for equilibrium 


air 


Other 

program variables are: 


DHRP 

P 2 Ah/Ap 


DP 

*i>/(e J 2 J 


DRP 

pM(K1/p>/3 P>2 


DRH 

p.vfoa/p.)>3h e ) 2 


ROR 

P„/Pj 
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E E (see eq. (A 10)) 

EMIN freestream Mach number, M 

A listing of subroutine SHOFF appears on the the following pages. . 
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Description of New Input Parameter 

i 

Only one additional input parameter is required, KE. The read state- 
ment, near the beginning of the main program, corresponding to format 
statement number 900 reads KE, where KE = 0 for normal-shock entropy and 
KE = 2 for variable entropy. These are the only allowable values for KE. . 

Description of New Output Parameters 

Near the beginning of the output from the program, the parameter KE, 
P2 (p 2 ), RH02 (p 2 ), H2 (h 2 ), B, BB (B) , RPER (R^, ) , RTB (R^ , DOR (A/R^,) , 

and NO . OF ITER (number of iterations for calculating the shock standoff 
distance) are printed. In addition, the following 8 parameters are printed 
for each position along a streamline: 

GBK local shock angle, f (5,3) 

GBL shock angle where streamline at boundary layer edge crossed the 

shock, ? (8,'!') 

PSH pressure aft of shock angle GBL, p g ^ (B ,40 (lb/ft^) 

PSHK pressure aft of local shock angle GBK, p gh (5,3) (lb/ft^) 

2 

PB surface pressure for shock solution only, p, (lb/ft ) 

D 

2 

PE surface pressure for calculating heating rates, p g (Ib/ft ) 

PSS(SHK) stream function at edge of shock layer, V ^ (slug/sec) 

PSS stream function at position where streamline at boundary layer 

edge crossed the shock, 4' (slug/sec) 

dL 
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TABLE 1. STAGNATION LINE STANDOFF DISTANCE COMPARED WITH MASLEN’S RESULTS 


p 2 /p » 

B 

A/R^ By Maslen's Method 
(Ref. 5) 

A/H t (Present Method) 

6 

2 

.06669 ± .001169 

. 06530 

20 

2 

.02462 ± .000328 

.02462 

50 

2 

.01098 + .0001195 

.01090 

200 

2 

.003046 ± .00002968 

.003022 

6 

1 

.1 ± 0.0 

.1010 

6 

2 

.06669 ± .001169 

.06530 

6 

3 

.05008 ± .001610 

.04795 

6 

5 

.03349 + .001963 

. 03104 

NOTE: 1) 

For all 

Calculations y =» 1.4, but the 

normal shock density 


ratio is assumed to be separately varied as shown. 


> 2) The B used herein is greater than Maslen's B; B = B^ ag ^ en + 



TABLE 2. DATA FOR BLUNTED 2:1 ELLIPTICAL CONE CASES 


P s (lb/in 2 ) 

T (°R) 
s 

M 

00 

o (lb/ft 2 ) 

CO 

T C°R) 
00 

1) 200 

1660 

10.04 

0.6395 

81.20 

2) 750 

1800 

10.13 

2.2766 

87.31 

3) 2000 

1900 

10.19 

5.9184 

92.00 



V (ft/sec) 

Re / ft 
00 

Re 

®,N 

Y 


1) 

4434.40 

.313 x 10 6 

2.61 x 10 4 

1.4 

.27354 

2) 

4639.42 

1.007 x 10 6 

8.39 x 10 4 

1.4 

.25010 

3) 

4790.61 

2.42 x 10 6 

2.02 x 10 5 

1.4 

.23470 
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TABLE 3. BODY SHAPE FOR SPACE SHUTTLE ORBITER WITHOUT CANOPY, 
LENGTH L “ 1.0 FT. 


Coefficients a^, for 


y, Windward 


Segment 

X i 

a 

o 

a l 

a 2 

a 3 

1 

.05 

.025 

.26795 

-4.0678 

30.1760 

2 

.10 

.032 

.08749 

-.4499 

.68725 

3 

.30 

.037 

-.010 

0 

. 

o 

4 

.50 

. 035 

-.010 

0 

0 

5 

.70 

■■1 

-.010 

0 

0 


y, Leeward 

Segment 

X 1 

a 

0 

a l 

a 2 

a 3 

1 

.05 

.025 

.26795 

0 

0 

2 

.10 

.0384 

.26795 

.12770 

-.88725 

3 

.30 

.090 

.21256 

-.22505 

-.4385 

1 

4 ; 

.50 

.120 

. 06993 

-.0366 

-.31525 

5 

.70 

.130 

.01746 

.02121 

-.07952 



z» 

Side 


Segment 

X 1 

a 

o 

a l 

a 2 

a 3 

1 

.05 

.025 

.26795 

-1.0446 

1.7120 

2 

.10 

.036 

.17633 

- .3383 

.90825 

3 

.30 

.065 

,150 

0 

0 

4 

.50 

.095 . 

.150 

- .02935 

.77175 

5 

.70 

. . 

.130 

.23078 

1.3697 

-1.7605 























































































Figure 2» Geometry of shock-wave angles r and a. 
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experiment, NASA TN D-5450 

variable entropy 1 streamlines from modified 
normal-shock entropy/ Ptessuraa 



Figure 5, Heat-transfer distribution on blunt 15° half-angle cone 
a ■ 20°, M *> 10.6, and Re „ = 3.75 x 10^. 
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Figure 9. Effect of Reynolds number on the heat-transfer distribution on 
blunted 2:1 elliptical cone, a = 30°, M = 10, 
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experiment, AEDC-TDR-64-19 
variable entropy 
normal-shock entropy 



Figure 11, Circumferential heat-transfer distribution on a blunted 
2:1 elliptical cone, a = 30°, M =10, Re .. ° 8.39 x 10 
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Figure 15. Heat-transfer distribution on typical space shuttle orbiter, M « 10, Re 
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Figure 17. . Circumferential heat-transfer distribution on typical space 

4 

shuttle orbiter, M = 10, a ° 30°, Re ,, 8 1.5 x 10 , 
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